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Determining a Prony Series for a Viscoelastic Material from Time 

Varying Strain Data 


ABSTRACT 

In this study a method of determining the coefficients in a Prony series representation 
of a viscoelastic modulus from rate dependent data is presented. Load versus time test data 
for a sequence of different rate loading segments is least- squares fitted to a Prony series 
hereditary integral model of the material tested. A nonlinear least squares regression 
algorithm is employed. The measured data includes ramp loading, relaxation, and unloading 
stress-strain data. The resulting Prony series, which captures strain rate loading and 
unloading effects, produces an excellent fit to the complex loading sequence. 


KEY WORDS: hereditary integral, viscoelasticity, weighted nonlinear regression, Prony 
series, multiple loading segments 


INTRODUCTION 


hi order to determine the time dependent stress - strain state in a linear viscoelastic 
material, under an arbitrary loading process, the deformation history must be considered. 

The time dependent constitutive equations of a solid viscoelastic material include these 
history effects. The load (stress) and displacement (strain) history, the loading rate 
(displacement rate), and time of load application on the specimen are all needed to determine 
the constants in the constitutive equations. A common form for these constitutive equations 


N 

employs a Prony series (i.e., a series of the form ^ (Xj 

i=l 


■e ~" r ' ). 


Creep and relaxation tests are most commonly used to determine the viscoelastic 
material properties, see Figure 1. In ideal relaxation and creep tests, the displacements or 
loads are applied to the specimen instantly. In the real test, especially for a large structural 
component, limitations of the testing equipment result in a relatively low strain rate and long 
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period of loading. The response during the period of loading is typically ignored, and only 
the data obtained during the period of constant displacement or constant load are used to 
determine the material properties . 1,2 Ignoring this long loading period and the strain rate 
effects in the data reduction introduces additional errors in the determination of the material 
properties. 

There are numerous methods for determining the Prony series from relaxation and/or 
creep data. An early method 3 involved constructing log-log plots of relaxation data in 
which straight line approximations for the data on the log-log graph yield the time constants 
(i.e. Tj ’s) from the slopes, and the exponential coefficients (i.e. a, ’s) are obtained from the 

intercepts. Other methods have also been employed. For example, Johnson and Quigley 4 
determined a relaxation time constant for a nonlinear model which is similar to a one-term 
Prony series model. They minimized a least-squares error measure, of the difference 
between the nonlinear model and measured data, by iteratively integrating (numerically) an 
internal variable equation. When attempting to determine relaxation time constants for 
higher fidelity nonlinear models, Johnson, et al. 5 employed trial and error procedures, 

Q 

similar to early linear methods, due to the complexity of the resulting nonlinear least- 
squares problem. More recently, a few authors 6 8 employed nonlinear optimization methods 
to obtain a high quality Prony series representation of relaxation data with a minimum 
number of terms in the series. The viscoelastic model can also be formulated in differential 
form. This is becoming popular recently 911 since the differential models can be effectively 
incorporated into finite element algorithms. When using these internal variable methods, 
each Prony series term is associated with a material internal state variable, hi the discrete 
finite element model, each term in the Prony series adds a substantial number of global 
variables. Thus, a short Prony series, which can accurately represent the material, is 
desirable. Nonlinear regression methods can help with determining a short and accurate 
Prony series. 

The purpose of this paper is to present a method for including the loading and 
unloading data, along with the relaxation data, in a nonlinear regression analysis to obtain the 
Prony series. The resulting viscoelastic material model is then capable of simulating the 
loading segments as well as the relaxation segments. This is an improvement when modeling 
hysteretic effects is important. The analytical solution for loading and/or unloading is 
determined herein and employed in a nonlinear regression analysis to determine the Prony 
series, hi addition, data weighting functions are investigated and are shown to improve the 
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fit in the beginning of relaxation period. Again, this method allows all the measured data to 
be included and results in an improved constitutive model. 


Hereditary Integrals for linear viscoelasticity 


1 2 

Detailed descriptions of linear viscoelasticity can be found in the literature. ’ An 
overview of linear viscoelasticity is provided here in order to introduce the hereditary 
integral method which is used below to determine the analytical solution for the loading and 
unloading segments. Linear viscoelastic constitutive models are represented by simple 
physical models composed of springs and dashpots. The spring is the linear-elastic 
component, and its constitutive equation is 

G — E • £ ( 1 ) 

The dashpot is the viscous component, and its constitutive equation is 


G —T\ ■ 


d£ 

St 


( 2 ) 


where r\ is the viscosity constant. Linear viscoelastic constitutive models are constructed by 
superimposing components with constitutive equations given by equations (1) and (2). 

Since the mechanical response of the dashpot is time dependent, the behavior of a 
viscoelastic material that is modeled by parallel and/or series combinations of springs and 
dashpots is also time dependent. 

The creep test consists of a constant stress, Oo , applied to a specimen for a period of 
time while its strain is recorded (Figure la). In a relaxation test, the specimen’s strain, £o , is 
held constant for a period of time while the stress is recorded (Figure lb), hi Figure 1, £o and 
Go are the initial strain and stress, respectively. For the relaxation test, a constitutive relation 
for the period of constant strain can be written as follows: 


< 7 ( 0 = 7(0 £ 0 ( 3 ) 

where Y(t) is a relaxation function. When the material is assumed to be a general Maxwell 
solid, the relaxation function is typically modeled with a Prony series as follows, 
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where: 


( 4 ) 


Y(t)=E 0 -0~lpi ■(\-e~ >lr ')) 

i = 1 


Pi is the i’th Prony constant ( i =12 ,. . .. ) 
h is the i’th Prony retardation time constant ( i = 1,. 2,.. .) 

E 0 is the instantaneous modulus of the material 
For time t = 0,Y(0) = E 0 and for t = , Y(oo) = E 0 ( 1-Ipt ). 

hi the case of a creep test, a creep compliance function, J(t), is defined as follows. 

£(t)=J(t)CF o (5) 

The compliance function is then determined by procedures analogous to those described 
above. 

To determine the stress state in a viscoelastic material at a given time, the 
deformation history must be considered. For linear viscoelastic materials, a superposition of 
hereditary integrals describes the time dependent response 1 . If a specimen is load free prior 
to the time t = 0, at which a stress, Oo + o(t), is applied the strain for time t > 0 can be 
represented as follows. 

where J(t) is the compliance function of the material and d(7(^) /d^is the stress rate. A 
similar equation can be used for the relaxation model to obtain the stress function introduced 
by an arbitrary strain function 1 

(7(l) = e l yY(t) + ^Y(t (7) 

where 7(t) is the relaxation function (Equation 4) and de ( 0 ld % is the strain rate. An 
example of applying hereditary integrals for a multiple loading segment process is shown in 
next section. 
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Hereditary integrals for a multiple loading process 


Hereditary integrals with Prony series kernels can be applied to model a loading 
process such as the one shown in Figure 2. The process in Figure 2 is divided into four 
segments for which strain and strain rate functions are defined. The functions are: 


e(t) = 


£ x t/(h -, o) 

t 0 <t<h 


e \Kh ~ f o) 


£ i 

h<t<t 2 

de 

0 

h < * < h 

-£yt /(t 3 - 1 2 ) 

h< t -h 

’ dt _< 

- 1 2 ) 

h <t -h 

0 

h < t - r 4 


0 

/j < t < < 4 


( 8 ) 


where £q = £ (0) = 0 and to = 0. 


For a material with a relaxation function in the form of a Prony series (Equation 4), 
the stress functions of the loading process can be derived as follows: 

Step 1. ( to < t< h) 

Substitute Equations 4 and first strain and strain rate functions of Equation 8 into Equation 7 
and obtain: 


<7(0 = YU) + 


n 




“0 + Jq£ 0 ‘(1“ 2 Pi •(!-<? 

/=i h 


_ Ep • £ \ 

h 

_ Eq -g i 


p£ + 'LPr*re 




t 


i 


Jo 

-t/ 




(9) 

where n is the number of terms in the Prony series. To simplify the expression, n will not be 
shown in following equations, pj and x, are the constants in the i-th term of the Prony series. 
Step 2. (ti<t< fe) 

Using the second strain rate function obtain: 
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< 7(0 = e w Y (0 + J ' -V (r - + \' h+ Y(t - #) ^ di f 
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= 0+El1± 
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Step 3. (t 2 <t< tj) 

The third strain rate function yields: 


CF(t) = e 0 -Y(t)+j , 0 '-Y(t-4)^^4+i^Y(l-4)^^d4 


de(g) 


d% 




+ji + r(*-f)^W 


= 0d 


'2+ ' ~ d% 

En • £ 


'0 C 1 


'1 


E 0 ■e l 


-(t- 


%-'LPiZ + 'LPi*i -e 
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+ 0 


Jo 


(h ~h ) 


%-'LPiZ + 'LPi* i -e 


-(t-a/" 1 ' 


Jt 2 


<t~h)y 


= ^ L (h-lp i h+lp i T i -e " ^-l Pl T r e~ A ) 


( 10 ) 


E 0 •£ 1 


-(t-hh 


(t~Y Pit + Y PFi ~h + Y Pit i - Y P/^- • e 7 7 ) 


(*3 ~ f 2 ) 

Note, tlie first portion of Equation 1 1 is equal to Equation 10. Thus, 


(11) 


~( t ~ t 2 ) 


<T(f) =<7 (0 - £ ° — (t - X Pit + X P/T; - t 2 + X P^2 - X PEi -e /Ti ) 
\h~t 2 ) 

where C 2 (t) is Equation 10 and it is function of time. 


( 12 ) 
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Step 4. ( f 3 < f < * 4 ) 

Similarly, the equation for the fourth step can be written as follows: 


cr(t) = CF 2 (t ) - 


Eg _£i 
(h ~h) 


(h-'LPih + 'LPfii 



~h 


+ 'LPih-'Lpfii 



) 


( 13 ) 


12 

A numerical example of a multiple loading segment process using MATHCAD 
software is shown in Appendix A. In the example, the stress function was calculated based 
on the strain and strain rate functions shown above, and it employed a two-term Prony series. 
The results of the viscoelastic analysis are shown in the stress-time and stress- strain plots. 
This worksheet can be used to generate data in a parametric study involving viscoelastic 
materials. Tire worksheet is also used as part of the weighted nonlinear regression algorithm 
as is shown in the following sections. 


Weighted Nonlinear Regression 

The Prony series coefficients and retardation times appearing in Equation 4 need to 

be determined in a regression analysis. Here, a standard nonlinear regression method (the 

13 14 

Marquardt-Levenberg Method ’ ) is used to perform the data fitting, hr the nonlinear 

2 

regression, an error function (y ) with respect to the unknown constants is defined as, 

9 N 

X (a) =t 

1=1 

where x, and yi are the experimental data, function y(xi;a) is the model to be fitted, and Oi is 

the standard deviation of measurement error of i-th data point. A set of unknown constants 

2 2 

(a) will be determined that minimize the error function % . Tire error function (y ) is 
approximated by its Taylor series with the quadratic form: 

X 2 (a)~c-d-a + y^-a-D-a (15) 

where c is a constant and d is the gradient of y~ with respect to the parameters a, which will 
be zero where y 1 is minimum. Matrix D is the second partial derivative of % 2 with respect to 
the parameters a. hritial trial values of a are specified and improved values are determined 


yi-y(Xiia) 


(14) 
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by the nonlinear regression algorithm. Iteration is continued until the error function, % , 
effectively stops decreasing. 

Since each Prony term includes two variables (p, and x,) and since the instantaneous 
modulus (Eq) must be determined, the total number of variables in the regression is 2/7+ 1 . 
Based on thermodynamic principles, several constraint conditions must be applied: 

P t > 0, I^<1, T,-> 0, E 0 > 0 (16) 

hi addition, the distribution of the standard deviation of measurement error (o) is not easily 
determined based on the error of data acquisition equipment and the error of test machine, the 
error is usually assumed to be uniform for all data points (o, = 1). As is well known, the 
viscoelastic effects are most significant at the beginning of the relaxation period, the fitting 
error in this region is significant. Since the percentage of the number of data points at the 
beginning of the relaxation period is less, the error function yj is dominated by a long 
uniform tail region of the relaxation period. To reduce the error and improve the fit at 
beginning of the relaxation period, a weight function (w = l/o, 0 < o <1) is used. The larger 
the weight factor at a data point, the better the curve fit the data point. There is no analytical 
method to determine the weight function, thus a trial- and-error method is used. The 
acceptance of the weight function is based on a graph of the data and regression model 
results. Weighted nonlinear regression requires more iterations than unweighted nonlinear 
regression, but it can provide a better fit to the experimental data in the region of most 
interest. 

Weighted Nonlinear Regression for Relaxation Test 

A three-point bending relaxation test of a composite material was performed 15 . The 
specimen (12 in. x 2 in. x 0.768 in.) was loaded in 22 seconds to a maximum deflection of 
0.103 in. at the middle of the span. Then the deflection was held for 11,711 seconds (Figure 

3). The Marquardt-Levenberg nonlinear regression method was applied by the commercial 

13 

software SigmaPlot to a hereditary integral model using two segments (loading and 
holding) to obtain the Prony series coefficients for the data. Since no analytical method 
exists to form the weight function, a trial-and-error method based on material properties was 
used to obtain a fit curve. The viscoelastic stress decays exponentially in the relaxation test, 
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therefore the viscoelastic effect is more significant during the loading period and at the 
beginning (< 30 seconds) of the holding period. The number of data points in these periods 

(<100) is much less than the number of data points in tail region of the relaxation period (> 

2 

5000). The error function (% ) will be dominated by a long uniform tail region of the 
relaxation period if a uniform weight function is applied. Therefore, a piecewise weight 
function was used to obtain better fits for these periods and improve the accuracy of the 
regression. Figure 4 shows the load relaxation at the beginning of the process. The dots 
represent the test data. Three regression results are shown. The dash-dot curve is the result 
of a regression analysis without the weight function (w/o WF) for a two-term Prony series. 
The long-dash curve is the result for a two-term Prony series with weight function number 1 
(WF1) shown below: 


150 <t <22 


w = < 


10022 < t < 50 
150 < t <11711 


(17) 


The short-dash curve is the result for a three-term Prony series with weight function 2 (WF2) 
as follows: 


W = < 


10t0<t <22 
10 6 /(7 2 )22 <t< 1000 
11000 <f <11711 


(18) 


Since an initial value is required for each of the variables, a trial data set based on the 
test data was assumed. The sum of the P constants should be about 0.09 since the load at the 
end of relaxation period is 9% lower than the load at the beginning. The retardation time 
constants can be set to powers of ten. As long as the initial trial values are reasonable, 
convergence will be achieved. 

As is well known, the regression data fits better in a particular region if the relative 
weight factor (weight factor / sum of weight factors) in that region is greater than the average 
value, hi Figure 4, the curves with weighted functions were closer to the data points near the 
beginning of the relaxation period. Since the weight factor at t = 22 seconds of function 2 
(=2066) is greater than the value of the weight function 1 (=100), the curve of function 2 is a 
better match to the data than curve of function 1 at the beginning of relaxation. 
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However, as shown in the Figure 3, weight function 2 does not fit the data as well as 
the other curves after 1000 seconds. It appears that function 2 was over weighted at the 
beginning of the relaxation and the relative weight factor at the other region was too small. 
The Prony constants of the regression are shown in Table 1. 


Table 1 . The Prony constants for the regression 



Modulus 

m 

Pi 

Tl 

(sec. ) 

P 2 

x 2 

(sec. ) 

P 3 

x 3 

(sec. ) 

w/o WF 

7259.4 

0.0259 

134.69 


4898.64 



WF 1 

7323.8 

0.0319 

57.23 


4040.36 



WF 2 

7450.9 

0.0262 

9.29 


80.65 

0.0460 

2028.03 


The results showed that a weight function should be selected based on the material 
properties and test data distribution. Properly selecting the weight function in the most 
significant region can improve accuracy of the regression. 

Weighted Nonlinear Regression for a Multiple Loading Process 

hi a recent study 15 , a thick composite panel responded viscoelasticly when it was 
tested, hi order to characterize the properties of the panel, a three-point bending test with 
multiple loading segment processes was performed on it. Since the stiffness of the panel was 
quite high, the test was conducted with a large hydraulic testing machine. High deformation 
rates were not available with this loading machine. The loading head displacement schedule 
of the machine is shown in Figure 5. The schedule is unlike a standard relaxation test. 
Though load relaxation periods exist, the time required to apply a full load was clearly long 
when compared to the relaxation periods. As mentioned in the previous sections, in order to 
include the data of loading and unloading segments, a nonlinear regression combined with 
the hereditary integral was used to simulate the data. 

Since the deformation was very small compared to the specimen dimensions a linear 
model was used for the calculations. The resulting displacements and the loads are linearly 
related to the strains and stresses. Thus, the displacement and load data were used as strain 
and stress data in the equations. 

The displacement (strain) schedule in Figure 5 was divided into 1 1 steps. The 
piecewise continuous function for the first 10 steps was assumed to be approximately linear, 
similar to the previous example (Equation 8). The data for step 1 1 was fit to a cubic 
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polynomial for a more precise regression model. The strain and strain rate functions were 
defined as follows, 


e(t) = 


-to) 


t 0 <t<t l 

h<t<t 2 


£; + G ' • t + b • t + • • • t^Q < t ^ t 




'll 


de 

dt 


e \! ( d\ ~ t o) t o <t -h 
0 




ci + 2 • b • t + • • • <t<t 


'n 


( 19 ) 


The parameters of the polynomial used in step 1 1 were s, = -6.47997, a = 0.022976, b 
= -2.516844e-5 and c = 8.772892e-9. 

The piecewise continuous stress function for first 10 steps, based on the hereditary 
integral, was derived by the same technique as shown from Equation 9 to Equation 13. The 
load function for step 1 1 was derived as follows: 

<r(r) = cr,o (0 + J ' o Y (t - 4) ^§^4 

= <7io(0 + J' £ 0 (l-ift(l-«' <K)/f, ))-(a + 2-&^ + 3-c^ 2 )-^ 

1=1 

= c7 io(0 + £o ‘[(<2 + b-£ 2 +c-£ 3 )- XPi + b-% 2 + c-£ 3 ) 

+ X Pre~ {t ~^ lTi • • (a + 2 • b • (% - Ti) + 3 • c • (<J 2 - 2 • ^ • r + 2 • r, 2 ))] ^ 

= (7 io(0 + Eo '[(l _ XP/)‘( a '( ,l_/l io) + ^'(^ 2 -tio) + c '( t3 _/ fo)) 

+ ^ Pi • Tj • (jd + 2 • b • (t — Tj ) + 3 • c • (J? — 2 • t • Tj +2 • Tj )) 

_ X/p • e -(r-rio)/r ' -Tj ■ (a + 2 ■ b ■ {t l0 - Tj) + 3-c-(t l0 2 -2-t l0 -z -\-2 -t 1 1 )) 

(20) 

where G, 0 (t) is the stress function from step 10. 

A piecewise weight function was generated in order to increase the accuracy of 
regression. Again, the weight function was determined by a trial-and-error method. After 
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several iterations, the constant weight factors used for each step were determine as: W(t) = 
[1.0, 5.0, 1.0, 1.5, 1.0, 7.0, 1.0, 1.5, 1.0, 5.0, 1.0] and the analysis results were shown in 
Figure 6. The weight factors were equal to 1.0 for the loading and unloading steps, and were 
greater than 1 .0 for the holding steps. 

The total number of variables in the Prony series material model is 2n+l, where n is 
the number of Prony terms. In this study, a two-term Prony series was sufficient to fit the 
data. The regression analysis results and test data are shown in Figures 6 and 7. Tire only 
difference between the results of the two analyses, shown in the figures, was the 
displacement function of step 1 1 . One of the analyses assumed that the displacement 
function is linear and other assumed it was cubic. Both model results agreed closely with the 
experimental data. Since the cubic polynomial described displacement in step 1 1 more 
precisely, the one with cubic displacement function fit better in step 1 1 than the one with 
lineal - function. The error function (% ) of the one with cubic function is 10% lower than the 
one with linear function. 


CONCLUDING REMARKS: 

A method of determining the coefficients in a Prony series representation of a 
viscoelastic modulus from rate dependent data has been presented. The hereditary integral 
method was employed to obtain an analytical representation of material response when it is 
subjected to rate dependent loading. The analytical representation was used in a nonlinear 
regression analysis, with measured data, to evaluate the Prony series constants. Several 
regression analyses were performed using different weight functions. For the data analyzed 
in this study, improved simulations of the hysteresis effects were obtained when the data at 
the beginning of each relaxation period was appropriately weighted. Note, the data analyzed 
here had loading and relaxation regions of similar length in time. Other weighting functions 
may be needed for different loading schedules. 

The method presented here provided a highly accurate representation of the material 
behavior in the rate dependent loading region. It can also represent the response of a 
viscoelastic material for other unique loading schedules. For example, it can be used for 
schedules in which the material is not allowed to relax between subsequent loading changes. 
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Appendix A Numerical Solution for Multiple Loading Segment Process 

The purpose of this appendix is to present an input file for MATHCAD which can be used to 
numerically compute results from Equations (9) to (13) for the multiple loading segment process. 
Two Prony terms are used in this simulation. The loading process is defined by Figure 2. hr the 
data file below, text with a bold font represents a comment and with a normal font represents a 
command. The result of the simulation is saved to a file (OUTPUT.PRN) which is shown at end 
of the appendix. 

Input file (Hereditary.MCD) 

The relaxation function of the material is as follows: 

Y(t) = E ( 1 - P1(1 - e - t/x1 ) - P2 (1-e - t/x2 )) 

Viscoelastic Material Constants: 

E:=10 9 PI := 0.2 xl := 10 P2:=0.1 x2 := 100 

where E is the modulus, P and x are Prony constants 
Loading time: (Seconds) 

delt := 5 holdt := 50 tO := 0 tl := tO + delt t2 := tl + holdt 

t3 := t2 + delt t4 := t3 + holdt 

Piecewise Strain Function: (Multiple Loading Process) t := tO . . t4 



0 otherwise t 


Hereditary Integral for the stress function 
Step 1 ( tO < t < tl) Loading to e = el 
t := tO.A 

o-l(r) :=E —-(t-Pl t + Pl tl-Pl Tl- exp(- — ) — P2 ■ t + P2 ■ t2 — P2 ■ t2 ■ exp(- — )) 
tl Tl t2 

Step 2 (tl < t < t2 ) Holding the load for 50 seconds 
t := fl..f2 

0-2(0 :=E-—(tl-Pltl + PlTl- exp(- t —^~) -Pl-tl- exp(- -)-P2-tl 

tl Tl Tl 

+ P2 ■ t 2 • exp(- - — — ) -P2-t2- exp( — — )) 
t2 f t2 
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Step 3 (t2<t<t3) Unload to e=e2 

t := t2..t3 


(72(t):=^^—^-(t-Pl-t + Pl-T\-P2-t + P2-T2-t2 + Pl-t2 + P2-t2 
t3-t2 

-Pl-tl- exp(- - — — ) - P2 ■ t2 ■ exp(- - — — )) 
rl t2 

ct 3(0 := <r2(r) + cr3_ l(r) 


Step 4 ( t3 < t < t4) Holding for 50s 
t := t3, t3 + 5..tA 

a A _ 1(0 := — ( - 2 — ■ ( t 3 -Plt3 + Plrl- exp(- ) - P2 • 1 3 + P2 • t 2 • exp(- - 12 

(63 - 12 ) rl t2 

+ PI • t2 + P2 ■ t2 - PI ■ Tl • exp(- ElE) -P 2 • t 2 • exp(- El^.)) 

rl t2 


(t 4(0:= a2(0 + (x4 _ 1(0 

Piecewise stress function 
t := t0..t4 


a(t) := 


(CTl(t)) if 0<t<tl 
(cr2(t)) if tl<t<t2 
(a3(t)) if t2<t<t3 
(c4(t)) if t3<t<t4 
0 otherwise 


Set up a matrix to print out of data 

X t,0 :=t V :=8(t) X t,2 =CT(t ) 


stress function 



Write the result file (OUTPUTPRN) 
WRITEPRN( output) := x 
(End of Hereditary. SV1CD) 


Output file (OUTPUT.PRN) 

Three columns of data (time, strain and stress) are 
included in the output file as shown below: 


(time) 

(strain) 

(stress) 

0 

0 

0 

1 

0.002 

1979653.653 

2 

0.004 

3921103.522 

3 

0.006 

5827816.446 

4 

0.008 

7702931.033 

5 

0.01 

9549288.871 

6 

0.01 

9389809.131 

7 

0.01 

9244678.865 


Stress vs. Strain 
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e(t) 


ie 


o(t) 

> 


Time 


characterization 
taxation test. 
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Strain 



Time (seconds) 

Figure 2. Multiple-segment loading process 
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Figure 3. The regression results of relaxation test 
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Load, lb. 
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Figure 4. The regression results at beginning 
of relaxation test. 


18 













0.00 0.05 0.10 0.15 0.20 0.25 0.30 


Displacement, in. 


Figure 7. Test and regression simulation of Normalized 
load-displacement for CAY panel 
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